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Critical Exponents for Supercooled Liquids 
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We compute critical exponents governing universal features of supercooled liquids through the 
effective theory of an overlap field. The static correlation length diverges with the Ising exponent; 
the dynamical correlation length grows more rapidly; and the relaxation time obeys a generalized 
Vogel-Fulcher-Tammann form. 



It is a longstanding problem to quantitatively charac- 
terize dramatic dynamical slowdown of supercooled liq- 
uids. Using the formalism developed in [![, we compute 
three critical exponents defined through 



^static 
^dynamic 

and log(r Q ) 



AT-", 
AT- 8 , 



(1) 
(2) 
(3) 



The results [cf. Eqs. (f2"2")) - (|24[) ] can be tested against ex- 
periments and simulations. 

Let us start by briefly reviewing the formalism. By 
doubly-rcplicating a system and defining a coarse-grained 
overlap field 4>, various observables of a single super- 
cooled liquid are imprinted onto the effective theory of 
an overlap field The static correlation length Static 
of the overlap field serves as a proxy for the typical size of 
cooperatively rearranging regions £crr required for the 
supercooled liquid to access distinct amorphous config- 
urations @]. It is then stipulated that relaxation of the 
overlap field constitutes a thermal activation process, me- 
diated by instantons. The size of this thermal instanton 
defines a dynamical length scale ^dynamic, a proxy for the 
typical size of mega rearrangements £m required for the 
system to erase coarse memories The energy cost 
of these instantons characterizes the metastable lifetime 
Tmetaj a proxy for the relaxation time scale r a . 

Envisioning that the number of accessible amorphous 
configurations becomes nonextensive near the Kauzmann 
temperature Q, we presume that the static correlation 
length diverges at a finite temperature To > 0. The ef- 
fective Hamiltonian near the critical temperature To is 
then argued to lie in the Ising universality class. Specifi- 
cally, as AT = T - To -> +0 and after sufficient coarse- 
graining, the effective Hamiltonian approaches that of 
the Ising model with its near critical temperature r ~ 

T i si „ g _ Rising AT < and its external magnetic 

field h ~ AT n (see FigHJ. Here n is required by analyt- 
icity to be a positive integer. An additional constraint on 
n comes from phenomenologically demanding that the ef- 
fective Hamiltonian hold a metastable state. As we shall 
see, this gives rise to the relation n > 2 in all spatial di- 
mensions d > 1. With the quick review passed us, we now 
determine the critical exponents for supercooled liquids. 




FIG. 1: A schematic picture of the line in the relevant param- 
eter space of the Ising universality class, onto which our effec- 
tive field theory projects. The dot denotes the unstable Ising 
fixed point where (r, h) — (0, 0) and the zigzag line marks the 
first-order phase transitions. The line approaching the fixed 
point is given by (r, h) ~ (-AT, AT") as AT -¥ +0 [J]. 



For d > 4, we have the effective Hamiltonian 
H[0(x)] 



k B T 



= / d d x 



(4) 



where u approaches a positive constant as AT — > +0. 
Before embarking on the instanton calculation, it is useful 
to rescale the variables as 



Irl 2 x and 6 = 



(5) 



so that the Hamiltonian takes the dimensionless form 



H[0(x) 
fc B T 

where 



= - I d d i 



and h = u 2 1 



(6) 



(7) 



The potential holds a metastable state for h < and 
this is guaranteed as AT — > +0 if and only if n > 2. 
Now, a thermal instanton 0i ns t (x) is a dominant saddle 
configuration which interpolates between the metastable 
minimum <^ me ta and the true minimum of the potential. 
It characterizes the leading metastable lifetime as Q 



^meta 



{AHA 

cxp I fcBT J 



with the activation energy 

AH a = H[0 inst (x)] - H[0 mo ta 



(8) 



(9) 
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As AT — > +0 with n > 2, the parameter /i approaches 
zero, making the calculation tractable as the thin-wall 
approximation [f| becomes asymptotically exact [l6| . 
In terms of the rescaled coordinates x, the instanton con- 
figuration has a thin-wall at radius 

i? inst ~ \h\~\ (10) 

Within the wall the field transitions between the two min- 
ima and outside is nearly constant. The activation energy 
of this configuration is given by 

AH^g-im-^-V. (11) 

Converting back to the original coordinates x, we thus 
obtain 

Static ~ ~ (AT)"* , (12) 



ensure the metastability condition by demanding n > 2. 
Near the critical point, we then obtain 

Static ~ h~ \r\-« ~ (ATy Ul , (19) 

Adynamic ~ ~ (AT)"^^ 2 "^^-^ , (20) 

and [13 

log(r Q ) ~ l^ffl-^-^ ~ (AT)-t (3 - A " W - 1)+(n - a)(,J - 1) > . 

(21) 

Finally, plugging in the known values of the Ising critical 
exponents [ll| yields 

v = vi (22) 
ps 0.63 for d = 3 
= 1 for d = 2, 



£dyn 



(AT) 



-{l + (n-2)} 



(13) 



and 



log (t q 



-(d-i) 



(AT) 



-{d-|+(n-2)(d-l)} 



(14) 

The previous mean-field scaling argument holds only 
for spatial dimensions d > 4. For d < 4, the would-be 
perturbative parameter g, governing higher order correc- 
tions, diverges near the critical point Q. We instead need 
to employ a renormalization-group scaling argument fl7| . 

For 1 < d < 4, let us again parametrize by r the rel- 
evant thermal direction and by h the relevant magnetic 
direction of the Ising universality class (cf. FigfT]). Under 
the renormalization group transformation, the thermal 
activation picture implies that 



I(r, h) = log (r mcta ) 



(15) 



has a simple scaling relation near the critical point. As 
we rescale the distance by a multiplicative factor 6, 



I(r,h) ps I{b 1/ui r,b^ /l,l h). 



(16) 



Here, l/v\ and A[ = /3j5i are the thermal exponent and 
the gap exponent of the Ising model, respectively. We 
rescale until b\ r = tmf with tmf fixed well away from 
the critical point, where the mean-field picture is valid. 
In particular, defining h e s = b^ h, 



(<*-!) 



J(r M F, keg) ~ \h e g 

Note that 

|^ ff |-M- A: l^l 



as 



leff 



(AT) r 



-a; 



(17) 



(18) 



approaches zero as AT — > +0 if and only if n > AJ: as 
Aj ps 1.56 for d = 3 and Aj = ^ for d — 2, we again 



V! + (2- AJ) + (n - 2) (23) 
1.07+ (n- 2) for rf = 3 

§ + (n - 2) for d = 2, 



and 



C = (2-Al)(d-l) + (»-2)(d-l) (24) 
ps 0.87 + 2(n-2) for d = 3 

= - + (n - 2) for d = 2. 
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The Ising exponent in the static correlation length has 
been observed in various numerical simulations fl2l . [l3j . 
In response to these, Langer proposed an Ising-like the- 
ory [14j. Work along this line may provide, not only 
evidence for the static part of our formalism, but also a 
complementary picture. 

As for the dynamics, we have 9 > is, indicating that the 
dynamic correlation length grows more rapidly than the 
static one. This "dccorrelation" of static and dynamic 
length scales is numerically observed in [TBI] . It would 
be interesting to see if our prediction for 9 can match 
experiments and simulations and, if it does, whether n is 
unique or labels different classes of supercooled liquids. 

The relaxation time obeys a generalized Vogel-Fulcher- 
Tammann form, 



exp 



fc B (T-T ) 



as T -> T , 



(25) 



with ( ^ 1. For d = 3, the minimal choice of n = 2 
yields £ ps 0.87, which is close enough to 1 that there is 
no apparent experimental conflict to our knowledge. It 
would be illuminating to check against extremely patient 
and accurate experiments in the future. 

Even if the formalism withstands rigorous experimen- 
tal tests, there will remain important questions. Do there 



3 



exist additional theoretical constraints on n? More fun- 
damentally, we assumed that the relaxation of the over- 
lap field proceeds through thermal activation. It would 
be exciting to figure out what kind of stochastic dynam- 
ics, Langevin or otherwise, governs the overlap field and 
whether it indeed gives rise to the thermal activation pic- 
ture. 
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